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ABSTRACT 

We report 350 /xm observations of 18 nearby luminous infrared galaxies 
(LIRGs), using the Submillimeter High Angular Resolution Camera II (SHARC- 
II) mounted on the Caltech Submillimeter Observatory (CSO) 10.4 m telescope. 
Combining our 350 /im flux measurements with the existing far-infrared (far- 
IR) and submillimeter (submm) data, we fit a single-temperature model to the 
spectral energy distributions (SEDs), and find the dust temperatures, emissiv- 
ity indices and far-IR luminosities having sample medians of Td = 39.4 ± 7.9 K, 
(3 = 1.6 ± 0.3 and LpiR = 10 ^^'^^"'^ Lq. An empirical inverse T^-jS correlation, 
best described by = [9.86 x 10^] is established for the local LIRG sam- 

ple, which we argue can be explained by the intrinsic interdependence between 
the dust temperature and grain emissivity index as physical parameters, as well 
as variations in grain properties in the the interstellar medium (ISM). 

Subject headings: dust — galaxies: infrared — submillimeter: galaxies 

1. Introduction 



Interstellar dust grains are small particles (~ 0.01 — 0.1 /im) that sparsely populate the 
ISM. Although interstellar dust accounts for only a very small fraction of the total mass in a 
galaxy, it plays critical roles in galaxy formation and evolution. Dust grains absorb strongly 
in the ultraviolet (UV) and optical, leading to a significant fraction of stellar radiation 
within a galaxy being absorbed. The warm grains subsequently emit strongly in the far- 
IR/submm, effectively down-converting the electromagnetic energy in various astrophysical 
environments. Observations of dust emission at far-IR and submm wavelengths are crucial 
in probing physical conditions and star formation activity. 

Emission and absorption line features are generally weak in the far-IR/submm (Blain 
et al. 2003), leading to smooth far-IR/submm SEDs dominated by dust thermal continuum 
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emission. However, modeling of the far-IR/submm SEDs observed in astronomical objects is 
far from trivial, as the observed emission spectrum is a complex function of radiative transfer 
as well as distributions in grain properties, such as composition, size and shape, that affect 
the way dust absorbs and emits radiation. In the simple case of a uniform grain population, 
dust thermal emission is well-approximated by a graybody (modified blackbody) function 
(Hildebrand 1983) 

= fiB.(Td)Q., (1) 

where Td is the dust temperature and Qi, is the absorption coefficient. Eq. ([T]) implicitly 
incorporates Kirchhoff's law, which ensures the equality between the emissivity and the 
absorption coefficient at all frequencies. In the far-IR/submm, = Qo (^)^, where Qo is 
the absorption coefficient normalized at some reference frequency and (3 is the emissivity 
index 0. There are more elaborate far-IR/submm SED models that attempt to account for 
multiple dust components and general optical depths. However, the applications of such 
complex models are impractical when the frequency sampling of the SED is limited. More 
importantly, Blain et al. (2003) show that the added complexity, even if feasible, generally 
does not lead to appreciable differences in constraining the observed far-IR/submm SEDs. 
Note that the SED parameters estimated by adopting the single-temperature SED model 
are effective, brightness-weighted average values from complex mixtures of dust grains with 
different temperatures, properties and optical depths in the ISM. 

Dust temperature profiles vary significantly for grain populations with distinctly differ- 
ent radii (i.e., a few submicrons versus tens of angstroms). The very small grains (a < 50 A) 
undergo large temperature fluctuations on absorption of a photon (Sellgren et al. 1985), 
while the larger "standard" grains reach equilibrium temperatures, as determined by ther- 
mal equilibrium between absorption of UV/optical photons and emission of far-IR/submm 
photons (Greenberg 1978), i.e., 

/ TT a^QijCUiydu = An / B,y{T^) na^ du , (2) 

JUV/optical Jfar— IR/submm 

where Uiy is the energy density of the incident radiation field. To derive the exact solution 
to Eq. ([2]), one would need to have accurate knowledge about the actual Qi, values at all 
frequencies - information that is lacking at present. At optical wavelengths, we know Q^, is 
relatively constant and near unity, and it is generally assumed that on average Q^, = 0.5 in 
the UV/optical in the numerical integration of the LHS of Eq. ([2]) (Martin 1978). While 
on the RHS of Eq. ([2]), Q,^ = Qo {u/uqY (by the definition of the emissivity index (3) at 



A blackbody has absorption coefficient Qj^ = 1 and emissivity index /3 = at all frequencies. 
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far-IR/submm wavelengths. For a uniform grain population, we thus have the approximate 
relation 

TdOc(-^)^, (3) 

where F is the integrated incident flux, defined as F = /uv/opticai '^^^^ characterizing the 
strength of the incident radiation field q 

Accurate knowledge of the emissivity index, which characterizes the frequency depen- 
dence of the absorption coefficient of interstellar dust, is crucial for determining important 
large-scale characteristics of dusty galaxies, such as dust temperature, dust mass, and the 
existence of multiple dust components. Given that the emissivity index is a complex function 
of various grain properties, its estimate also offers insights into grain properties, particularly 
grain structure (Seki & Yamamoto 1980; Yang 2006) and formation of ice mantles (Aannes- 
tad 1975). While theory predicts (3 values of 1 or 2 for interstellar dust in the long- wavelength 
limit, observations of astronomical objects in the far-IR/submm have revealed a wide range 
of (3 values. Lower emissivity indices have been observed in active galactic environments, 
such as circumstellar disks and warm molecular clouds, as well as large-scale dust distribu- 
tions in infrared luminous dusty galaxies (Knapp et al. 1993; Blake et al. 1996; Dunne et 
al. 2000). Far-IR emissivity that has a frequency dependence steeper than quadratic has 
also been observed in some galactic sources (Schwartz 1982; Lis et al. 1998). In particular, 
Dupac et al. (2003) reported a wide range in /3 values, as well as a negative temperature 
dependence of /3, as observed in a large sample of Galactic molecular clouds. 



2. Sample and Observations 

Sources in this sample are mostly selected from the LIRG sample presented by Wynn- 
William & Becklin (WWB, 1993). The WWB sample, not being a strictly defined complete 
sample, includes 19 sources that were selected from the IRAS Point Source Catalog for their 
high luminosities and high fiux densities at 12 and 25 /im. 16 sources here are selected from 
the WWB sample by requiring Sioo^tm > 20 Jy. In addition, two ULIRGs - UGC5101 and 
Mrk231 - are added to the list to better understand nearby infrared galaxies at the high 
luminosity end. In total, our sample includes 18 nearby LIRGs in the redshift range of 
0.003 < z < 0.042. This sample, with its selection based on high luminosity and high fiux 
densities at the IRAS wavebands, is likely to have an enhanced fraction of galaxies that are 



2 p ^ 8^ (kTd/hz.o)"+'' r(4 + P) C(4 + /3), where r(z) = /~ t^-^e"' dt and C(s) = /o°° dt 
are the Gamma and Riemann C, functions respectively. This formula is equivalent to that given by De Breuck 
et al. (2003). 
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more active and contain warmer dust as compared to the bulk of the submm population 
at similar redshifts. A wealth of photometry data at far-IR/submm/mm wavelengths is 
available from the literature for a majority of the sources, which enables us to estimate both 
dust temperature and emissivity index in the SED fitting. Table [T] and Table [2] list the 
sample and the existing far-IR/submm/mm photometric data. 

The observations were made during a series of SHARC-II (Dowell et al. 2003) observing 
runs from January 2003 to September 2004 at the CSO under moderate weather conditions 
(0.07 < T225GHZ < 0.10). Integration time varied from source to source, depending on the 
brightness of each source at 350 /im, and the atmospheric opacity and sky variability. On 
average, 0.5 hours of integration time was spent on each source. All measurements were made 
by scanning the bolometer array in a Lissajous pattern centered on each source; scans were 
taken without chopping the secondary mirror. Pointing and calibration scans were taken 
on an hourly basis on strong submillimeter sources. The absolute calibration was found to 
be accurate within 20% error margin. Raw data were reduced using the software package 
CRUSH (Kovacs 2006), to obtain a 350 /im map and flux density for each source. All of 
the 18 sources were clearly detected at levels of S/N > 10, and Fig. [T] shows the SHARC- 
II 350 /zm maps clipped to have at least half of the maximum exposure. Table [3] lists the 
350 /im flux densities and flux errors. 



3. Spectral Fits and Derived Quantities 

We adopt the single-temperature, optically-thin far-IR/submm SED model, as given 
by Eq. ([T]). A flux error at a uniform 20% is assumed at each data point, hence each flux 
measurement is weighted equally in the nonlinear least squares fit. All but four galaxies 
(NGC2388, NGC4194, NGC4818 and NGC5135) in the sample have photometry data at 
far-IR/summ/mm wavelengths in addition to 60, 100 and 350 /im, thus the dust temperature 
and emissivity index can both be treated as free parameters. However, the SED fitting proce- 
dure introduces a significant degeneracy between the estimated values of Tj and (3 (Blain et 
al. 2003) - a problem that persists even when the frequency sampling of the far-IR/submm 
SED is highly redundant. To achieve redundancy in the cases of NGC2388, NGC4194, 
NGC4818 and NGC5135, where the sampling of the SED is limited to three wavelengths, 
we assume /? = 1.5 - a value most commonly assumed and observed. 

The best-fitting dust temperature has a sample median and standard deviation of Td = 
39.4 ± 7.9 K and ranges within 29.2 < Td < 62.9 K, with two galaxies (NGC4418 and 
Mrk231) showing Td > 50 K. The emissivity index ranges within 0.9 < (3 < 2.4, and has a 
sample median of (3 = 1.6 ± 0.3. The far-IR/submm SED fits of our local LIRG sample are 
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shown in Fig. [2l the distributions of and (3 are shown in Fig. [31 

The global far-IR luminosity (Lpu^) and dust mass (M^) are calculated from 

/•lOOO Atm 

LpiR = AtiDI / S^du (4) 

J 40 fira 

Ma = (5) 

where Dl is the luminosity distance, and k^, = |^ is the dust mass absorption coefficient, 
varying as oc u^. We adopt dust parameters given by Hildebrand (1983): Qi25/im = 7.5 xlO""^, 
p = 3.0g/cm^ and a = 0.1 /im, which implies Hi25fim = 1.875 (kg/m^)~^. While there is 
significant uncertainty associated with k,, due to the lack of accurate knowledge of interstellar 
dust properties, Lpi^ is tightly constrained when a satisfactory SED fitting is achieved. The 
LpiR and Md values inferred for our local LIRG sample range within 10®'^ < Lpi^ < 10^^'^ 
L© and 10 < Md < 10 M©, and have sample medians of LpH, = lO^ ^iO-^L© and 
Md = 10 ''•^^^■^ Mq. Thus, our sample spans more than two orders of magnitude in both 
far-IR luminosity and dust mass. Table [3] lists the fitted values of Td and f3, along with the 
LpiR and Md estimates for our local LIRG sample. 



4. Td-/5 Relation 

4.1. Correlation Coefficient &; Simulations 

A negative correlation between the dust temperature and emissivity index clearly emerges 
from the T^-jS scatter plot for our local LIRG sample (Fig. H]). Note that, in the investiga- 
tion of the Td-/5 relation, we only include the 14 sources for which the frequency sampling is 
sufficient (> 4) to allow simultaneous estimates of Td and (3. The non-parametric correlation 
coefficient d between the dust temperature and emissivity index estimated for our 14 local 
LIRGs is calculated to have a value of p^p = —0.79, and it is non-zero at the significance 
level of (1 — p), where p ~ 7.43 x 10"'^ (Yang 2006). We thus conclude the inverse Td-/3 
correlation observed in our local LIRG sample is statistically significant. 

However, we note that fitted values of Td and (3 are simultaneously obtained from the 
SED fitting, in which a strong and negative Td-/3 degeneracy is present and could potentially 



•^Without prior knowledge of the true distributions of dust temperature and emissivity index, the non- 
parametric (rank) correlation coefficient between these two variables provides a more robust correlation test 
than the parametric correlation coefficient (Press 1999). 
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introduce an artificial inverse Td-/3 correlation. We address this problem by repeatedly 
performing SED fittings on simulated photometric data. Our simulation method, described 
as follows, is a modification and extension of the simulation procedure performed by Dupac 
et al. (2002). First, we are primarily concerned with the non-parametric Tj-/? correlation 
coefficient; second, the effects of the Td-P degeneracy on the observed Td-/? correlation can 
be accurately quantified. 

Two sets of independent, uniformly distributed random numbers are generated within 
the ranges of 20 — 60 and 1.0 — 2.5, and they are taken to be the simulated values of Td and 
f3 respectively for each of the 14 local LIRGs. Furthermore, we assume flux errors to be 
normally distributed and at a constant level of 20% at all wavelengths. Using the simulated 
values of Td, (3 and measurement errors, simulated photometric data are calculated at all far- 
IR/submm wavelengths where observed data are available for each of the 14 sources in our 
local LIRG sample. SED fitting of the simulated photometric data for each source follows, 
yielding fitted values of Td and f3. We calculate non-parametric correlation coefficients 
between the "real" (simulated) as well as the fitted values of Td and (3 in the sample. 

The above steps are repeated many times (~ 10000), and the probability distribu- 
tion function of the non-parametric correlation coefficient between the dust temperature 
and emissivity index is thereby established. As shown in Fig. O pnp of the simulated Td- 
j3 show a distribution well approximated by a Gaussian function ~ N(— 0.01, 0.29), as can 
be expected from the central limit theorem given that the simulated Td and (3 are inde- 
pendent random variables. By contrast, pnp of the fitted Td-/3 roughly follows a Gaussian 
distribution ~ N(— 0.14, 0.29), suggesting the SED fitting procedure typically introduces a 
non-parametric Td-/3 correlation coefficient at the level of pnp ~ —0.14 when Td and (3 are 
in fact uncorrelated, given the SED frequency samplings of our 14 local LIRGs in the far- 
IR/submm/mm. Under the normal distribution ~ N(— 0.14, 0.29), the observed Td-/? corre- 
lation coefficient, pnp = —0.79, is significant at the level of (1 —p*), where p* ~ 1.20 x 10~^. 
Hence, we deem the negative Td-/3 degeneracy in the SED fitting insufficient to explain the 
inverse Td-/3 correlation observed in the local LIRG sample and conclude that the observed 
inverse T^-jS correlation is real. 

4.2. Functional Fitting 

Based on physical considerations, particularly Eq. ([3]), we fit the inverse Td- f3 correla- 
tion observed in our local LIRG sample using the function 



= Ci '^z+zs 



(6) 
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The best fit to Eq. ([6]) is 



Td 



[9.86 X 10^] , 



(7) 



which provides a good fit to the dust temperatures and emissivity indices estimated for our 
14 local LIRGs, as shown in Fig. HI We also try two alternative functional forms in which 
the dust emissivity index is a power or hyperbolic function of the dust temperature, i.e. 



The best fits to Eqs. & and (ED are /? = [82.37] Td and f3 = [62.89]/Td respectively, 
and yield Td-P curves nearly identical to Eq. Thus, we consider Eq. ([71) as the best 
description of the Td-/3 relation observed in the local LIRG sample. This Td-/3 relation is 
also tested in the SED fittings of a ULIRG sample at intermediate redshifts (0.1 < z <, 1.0, 
Yang 2006), for which the incorporation of the Td-/3 relation (Eq. [7j) leads to reasonable 
SED fits for the vast majority of the sources. This provides strong support for an extension 
of the Td -jS relation, as derived for our local LIRG sample, to dusty galaxies in the more 
distant Universe. 

For comparative purpose, we also plot the Td-P relation given by Dupac et al. (2003), 
as derived from multiband submm observations of a large sample of molecular clouds in 
the Galaxy (Fig. H]). These authors reported a wide range of Td and (3 values estimated 
from SED fittings and found an inverse Td-/3 correlation best fitted by a hyperbolic func- 
tion P = 4+0^08 Td • Clearly, the Td-/3 correlation seen in the Galactic molecular clouds 
is significantly different from, in fact much flatter than, that observed in our local LIRG 
sample, despite the agreement in qualitative trend. We note that the Td and /3 estimates 
for the molecular clouds are sometimes given by SED flttings without flux measurements at 
wavelengths shortward of the SED peaks, which would otherwise have greatly improved the 
accuracy of the Td estimates. The difference in the shape of the observed Td-P correlations, 
as derived for astronomical samples of differing nature, scales and environments, is intrigu- 
ing. We argue that variations in the Td-/3 relation are not unreasonable, given that various 
physical characteristics, such as the radiation fleld, grain properties and grain distributions, 
could all impact the large-scale dust temperature and emissivity index, and thus the way 
in which Td and /3 relate to each other. There is considerable observational challenge cur- 
rently facing us in the task of achieving sufficient far-IR/submm/mm SED samplings, for a 
large number of dusty objects over a wide range in redshift, physical scale and environment. 
Dedicated observational efforts in the future would allow us to more precisely quantify the 
Td-/3 relation, discover its full range and investigate its evolution. 



ciTd' 
1 



(9) 



8 



ci + C2 Td ■ 
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4.3. Laboratory Results 

Laboratory measurements of the absorption coefficients of interstellar grain analogs at 
submm/mm wavelengths have revealed an intrinsic temperature dependence of grain absorp- 
tion coefficient and emissivity index (Agladze et al. 1996; Mennella et al. 1998). Agladze et 
al. (1996) measured absorption spectra of crystalline and amorphous grains at wavelengths 
ranging between 700 fim and 2.9 mm at temperatures between 1.2 and 30 K. The emissivity 
index was found to range between 1.2 and 2.5, and show a negative temperature depen- 
dence for amorphous grains at temperatures above 10 K. Mennella et al. (1998) measured 
absorption coefficients of crystalline and amorphous dust analog grains between 20 fim and 
2 mm over the temperature range 24—295 K. The temperature and wavelength domains 
covered in this set of laboratory experiments are probably more relevant for understanding 
astronomical observations at far-IR/submm wavelengths. The emissivity index, as well as 
its temperature dependence, was found to be sensitive to grain material and grain struc- 
ture, and variations in the absorption coefficient and emissivity index were found to be more 
prominent in amorphous grains than in crystalline grains. These authors found a monotonic 
increase of the absorption coefficient with dust temperature, while the emissivity index was 
found to range from 0.5 to 2.3 and decrease as the dust temperature increases. 

4.4. Explanations 

We conclude the inverse Td-P correlation observed in our local LIRG sample is statis- 
tically significant and physically real, and offer three possible explanations. 

First, Eq. ([3]) tells us that the equilibrium dust temperature displays a negative depen- 
dence on the emissivity index, for a given radiation field and grain material. 

Second, the grain emissivity index has intrinsic negative temperature dependence, as 
evidenced by laboratory measurements (§ 14. 3p . This effect can be understood by considering 
the Schlomann (1964) model which treats an amorphous solid as an ionic crystal with a ran- 
dom charge distribution. Within this framework, the absorption spectrum is dominated by 
photon-to-phonon coupling at all frequencies. We argue that the observed negative temper- 
ature dependence of grain emissivity index is attributed to the decreasing grain dimension 
(D) at higher temperatures (i.e., grain structure becomes more "open") 0, which results in a 



"^Here we adopt the concept of fractal dimension (Hausdorff dimension), defined as D= logp/pN, where 
P/p is the size ratio of self-similar fractal objects and N is the number of smaller units that can fit into the 
larger unit. Structures with lower fractal dimension are more "open". 
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flatter frequency dependence of the phonon mode density and consequently lower emissivity 
indices at higher dust temperatures (Yang 2006). 

Third, an inverse T^-fi relation can be caused by mixtures of non-uniform grain pop- 
ulations, particularly those involving very small grains or grains covered with ice mantles. 
/3 values as high as 3.5 have been predicted for dust grains covered with ice mantles (Aannes- 
tad 1975), whereas very small grains are expected to have /3 ~ 1 (Seki & Yamamoto 1980). 
At the same time, low dust temperatures favor formation of ice mantles and lead to strong 
emission at longer wavelengths, while high dust temperatures (or temperature fluctuations) 
in very small grains prohibit accretion of ice mantles and give rise to strong emission at 
shorter wavelengths. As a result, the dust temperature and emissivity index observed in the 
large-scale tend to move in the direction of low-Td-high-/? when the overall dust emission is 
dominated by that of cold dust possibly covered with ice mantles, while Ta and /3 shift in 
the direction of high-Td-low-/3 when very small grains dominate the overall emission. 

4.5. Applications 

A Td-P relation, if established over a wide range of dusty galaxies, would be extremely 
useful for modeling their observed far-IR/submm SEDs, as it effectively reduces the number 
of variables (therefore the number of required flux measurements) by one in the SED fltting. 
This effect is especially desirable for studying faint, distant submm sources and interpreting 
results from single submm wavelength imaging and deep fleld surveys, given the considerable 
observational challenge typically faced by deep submm observations. For a given photometric 
dataset, a T^-/? relation either allows the estimate of one more parameter in the SED fltting, 
or improves the uncertainties associated with the fltted variables, as exemplifled by the SED 
modeling of a sample of intermediate-redshift ULIRGs (Yang 2006), for which the sampling 
of the far-IR/submm SED is hmited to three wavelengths. By using the T^-P relation, 
both the dust temperature and emissivity index can be fltted, which represents a major 
advance over the more conventional and commonly utilized SED modeling method, which, 
without the application of a Td-/3 relation, would only flt for one of the two very important 
parameters. 

5. Conclusion 

We have observed and detected 18 local LIRGs (0.003 < z < 0.042) at 350 //m, using the 
SHARC-II camera at the CSO. The acquired 350 iim data, in combination with pre-existing 
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far-IR/submm/mm photometric data, lead to accurate estimates of SED properties of these 
galaxies. We find sample medians of dust temperature (T^ = 39.4 ± 7.9 K), emissivity index 
{p = 1.6 ± 0.3), far-IR luminosity (Lp„, = 10 ^^-^^o-e Lq) and dust mass (Md = lO^-^^o-^ Mq). 
An inverse Tj-/^ correlation, best described by = [9.86 x 10^] 4:^3+^, is established for our 
local LIRG sample. This effect is most likely caused by intrinsic interdependence between 
the dust temperature and emissivity index as physical parameters, as revealed by laboratory 
measurements (§ 14. 3p . as as well as variations in grain properties in the ISM. The estab- 
lishment of a Td-(3 relation is extremely useful in the far-IR/submm SED fittings of dusty 
galaxies, as it serves to effectively reduce the number of free parameters by one in the SED 
function. 

We would like to thank J. Zmuidzinas, A. Blain, T. Soifer, T. Greve and C. Dowell for 
helpful discussions. We are grateful to the referee for improving this paper through thoughtful 
comments. The CSO is supported by the NSF fund under contract AST 02-29008. 
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Table 1. The local LIRG sample selected for SHARC-II observations 



Source Name 


RA 




Dec 




z 


Seo 


Sioo 




J2000 




J2000 






Jy 


Jy 


NGC520 


011i24m35. 


Is 


+03d47m33s 


0, 


.007609 


30, 


.37 


46.15 


NGC1614 


041i33m59. 


,8s 


-08d34m44s 


0, 


.015938 


33. 


.02 


34.35 


NGC 2339 


071i08m20. 


,5s 


+18d46m49s 


0, 


.007358 


18, 


.45 


31.46 


NGC2388 


07h28m53, 


,4s 


+33d49m09s 


0, 


.013790 


17, 


.03 


25.33 


NGC 2623 


08h38m24. 


,1s 


+25d45ml7s 


0, 


.018463 


23, 


.80 


26.66 


UGC5101 


09h35m51. 


,6s 


+61d21mlls 


0, 


.039390 


12, 


.24 


20.25 


NGC 4102 


12h06m23. 


,1s 


+52d42m39s 


0, 


.002823 


47, 


.04 


73.84 


NGC 4194 


12hl4m09. 


,6s 


+54d31m36s 


0, 


.008359 


24, 


.09 


26.06 


NGC 4418 


121i26m54. 


,6s 


-00d52m39s 


0, 


.007268 


45. 


.58 


31.99 


Mrk231 


121i56ml4. 


,2s 


+56d52m25s 


0, 


.042170 


33 


.28 


30.29 


NGC 4818 


121i56m48. 


,9s 


-08d31m31s 


0, 


.003552 


20. 


.32 


26.72 


NGC 5135 


13h25m44. 


,0s 


-29d50m01s 


0, 


.013716 


16. 


.93 


30.72 


Mrk273 


131i44m42. 


,1s 


+55d53ml3s 


0, 


.037780 


23 


.09 


21.97 


NGC 6000 


15h49m49. 


,5s 


-29d23ml3s 


0, 


.007315 


36, 


,48 


51.76 


NGC 6240 


16h52m58. 


,9s 


+02d24m03s 


0, 


.024480 


22, 


.87 


27.21 


IC5135 


211i48ml9. 


,5s 


-34d57m05s 


0, 


.016151 


16, 


.36 


26.19 


NGC 7469 


231i03ml5. 


,6s 


+08d52m26s 


0, 


.016317 


26, 


.49 


36.15 


Mrk331 


231i51m26. 


,8s 


+20d35ml0s 


0, 


.018483 


18, 


.30 


23.44 



Note. — 



Given by NASA Extragalactic Database (NED). 



Table 2. Other far-IR/submm/mm photometric data available for the local LIRG sample 



Source Name 


far-IR/submm/mm Photometric Data 




Jy 



NGC520 
NGC1614 
NGC2339 
NGC2388 
NGC 2623 
UGC5101 
NGC 4102 
NGC 4194 
NGC 4418 

Mrk231 



NGC 4818 
NGC 5135 
Mrk273 

NGC 6000 
NGC 6240 

IC5135 
NGC 7469 
Mrk331 



SgSOum 
SgSOum 
S450um 

SsSOum 

S450um 

c(b) 
'-'450um 

c(b) 
'^450um 
q(b) 
'-'1.3mm 

Sl20um 
S450um 

q(b) 

1.3mm 



Sl20um 

S450um 

c(b) 
'-'450um 

Sl20um 

S450um 

Sl.3mm 

SsSOum 

SsSOum 



0.325 ± 0.050, sSi^^, = 0.050 ± 0.025 
0.219 ± 0.032, Sg^^ - 0.011 ± 0.003 
0.849 ± 0.078, Ssooum = 0.079 ± 0.015, Si.amm 



0.019 ±0.002 



0.091 ±0.014 

1.43 ± 0.304, Sgooum = 0.143 ± 0.025, Si.i^^ = 0.068 ± 0.012 



1.79 ±0.450, s£Lm 

j(b) 



0.212 ± 0.028, SS^^ = 0.055 ± 0.007, Si";^^^ = 0.050 ± 0.006 



(b) 



0.240 ± 0.020, Sgsoum = 0.255 ± 0.037, Si' 



(b) 



1mm 



0.085 ± 0.009 



1.34 ±0.350, S^"ooum 
0.072 ±0.012 

24.32 ± 30%, Sisoum = 14.74 ± 30%, Sigoum = 9.75 ± 30%, Saooum = 6.88 ± 30% 

0.513 ± 0.134, S^SIj^,^ = 0.085 ± 0.012, sf}^^ = 0.045 ± 0.012, Si.25mm = 0.029 ± 0.008 

0.038 ± 0.008 



= 19.99 ± 30%, Sisoum = 13.10 ± 30%, Sigoum = 8.69 ± 30%, S2ooum = 7.40 ± 30% 

= 0.707 ± 0.200, Ssooum = 0.084 ± 0.022, Si.ii„m = 0.051 ± 0.013, S^^'l^^ = 0.020 ± 0.002 

= 1.460 ± 0.650, SS^^ = 0.115 ± 0.027, sj^}^^ = 0.055 ± 0.008, Sf]^^ = 0.018 ± 0.008 

= 25.90 ± 30%, Sisoum = 18.91 ± 30%, Sigoum = 12.73 ± 30%, Saooum = 9.00 ± 30% 

= 1.00 ± 30%, Sssoum = 0.150 ± 30%, Sg^^ = 0.021 ± 0.010 

= 0.033 ±0.004 

= 0.192 ±0.027 

= 0.132 ±0.025 



Note. — Given by NED unless noted otherwise. 

References. — (a) Benford 1999; (b) Roche & Chandler 1993. 
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Tabic 3. 350 //m flux and derived properties for the local LIRG sample 



Source Name 


S350 


C350 


Td 


^(Td) 


/5 


a{(3) 


log LpiR 


logMd 




Jy 


Jy 


K 


K 






L© 


Mq 


IN Vjr OZi\J 


O.Ul 


n 09 


/ .0 




i-.O 


9 
u.z 


1 n 7fi 


fi QA 


NCIC 1 fil ZL 


1. 10 


n 


An 7 


n 


1 Q 


n 9 


1 1 '^7 
1 1 .0 1 


7 '?n 


IN Of v_/ ZOOy 


1. ( 


n OA 


^^9 1 


9 


9 1 


n 9 




7 ni 




L.O 1 


n 1 n 


^^7 Q 


9 
o.z 


i.O 




i i .uo 


7 91 






u.uo 


41 8 


u.u 


l.U 


u.o 


1 1 ^^8 
1 1.00 


7 ^^9 


UGC5101 


2.20 


0.10 


34.1 


3.6 


1.6 


0.2 


11.86 


8.27 


NGC4102 


4.38 


0.04 


32.1 


3.0 


2.1 


0.2 


10.10 


6.56 


NGC4194 


1.19 


0.03 


45.3 


4.7 


1.5 




10.68 


6.47 


NGG4418 


2.03 


0.07 


62.9 


12.4 


0.9 


0.2 


10.81 


6.14 


Mrk231 


1.73 


0.07 


51.3 


5.6 


1.4 


0.1 


12.24 


7.79 


NGC4818 


2.73 


0.05 


36.0 


3.2 


1.5 




9.91 


6.20 


NGC5135 


3.80 


0.07 


33.0 


2.6 


1.5 




11.10 


7.59 


Mrk273 


1.34 


0.08 


42.7 


4.0 


1.6 


0.1 


11.98 


7.88 


NGC6000 


4.26 


0.20 


29.2 


2.6 


2.4 


0.2 


10.81 


7.48 


NGC 6240 


1.10 


0.05 


39.4 


3.4 


1.7 


0.1 


11.64 


7.69 


IG5135 


2.16 


0.10 


35.5 


3.9 


1.6 


0.2 


11.17 


7.49 


NGG7469 


2.23 


0.20 


42.3 


6.9 


1.3 


0.3 


11.35 


7.31 


Mrk331 


1.69 


0.04 


41.5 


6.5 


1.4 


0.3 


11.29 


7.30 



Note. — (3 =1.5 is assumed in the SED fitting for NGC 2388, NGC4194, 
NGC 4818 and NGC 5135. The luminosity distance used in the calculation of 
the far-IR luminosity is derived from the redshift, by assuming a flat cosmology 
with Ho = 70kms-i Mpc'^ = 0.30 and Qa = 0.70. 



- 16 




1100 50 -50 -100 1100 50 , -50 -100 MOO 50 -50 -100 MOO 50 -50 -100 
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Fig. 1. — SHARC-II 350 /im maps of galaxies in the local LIRG sample. In each map, the 
green plane represents the part of bolometer array which has at least half of the maximum 
exposure, and the red peak in the middle of the array indicates a detected source (S/N > 10). 
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Fig. 2. — SED fits for galaxies in tlie local LIRG sample. 
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Fig. 2. — Continued 
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Fig. 2. — Continued 



A (um) 



A (um) 



A (um) 



1.0E + 02 



1.0E+01 



1.0E+00 



1.0E-01 



l.OE-02 



1000.0 



NGC6000 
T= 29. 8K 
13= 2.3 

100.0 




1.0E+01 



l.OE+00 



1.0E-01 



l.OE-02 



A (um^ 



1000.0 



IC5135 
T= 37. IK 
1.5 

100.0 




1.0E+01 



l.OE+00 



1.0E-01 



l.OE-02 



A (um) 



MRK331 
T= 39. 8K 
|8= 1.4 



1000.0 



100.0 



A (um) 



1000.0 



1000.0 



1000.0 



-20- 




Fig. 3.- 



Histograms of Td and (3 derived for the local LIRG sample. 
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Fig. 4. — The inverse Td-/5 correlation observed in our local LIRG sample when the far- 
IR/submm/mm SEDs are sampled at no less than four wavelengths. The dashed line rep- 
resents the best fit to the T^ and (3 values derived for the 14 local LIRGs. For comparison, 
the Td-/3 correlation derived for the Galactic molecular cloud sample (Dupac et al. 2003) is 
shown by the dashed-dotted line. 
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Fig. 5. — a): The distribution of p^p calculated for the uncorrelated, simulated Td-f3 is well 
approximated by a Gaussian function ~ N(— 0.01, 0.29), as shown by the solid line, b): 
The distribution of pnp calculated for the fitted Td-/9, from SED fittings of the simulated 
photometric data, is well approximated by a Gaussian function ~ N(— 0.14, 0.29), as shown 
by the dashed line. Under this distribution, the pnp of the fitted Tj-/?, from SED fittings of 
the observed photometric data, is clearly significant at a value of pnp = —0.79 (marked by 
the sticker) . 



